Self-heating in small mesa structures 
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We study analytically and numerically a problem of self-heating in small mesa structures. Our 
results show that the self-heating is proportional to a characteristic in-plane size of the mesa. Ex- 
perimental data for small high-T c superconductor Bi2212 mesas are in qualitative agreement with 
our calculations. We estimate the self-heating in Bi2212 mesas with different sizes and demonstrate 
that the self-heating can effectively be obviated in small mesa structures. 
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Progress in microelectronics results in a continuous de- 
crease of sizes of individual elements. Often, an active 
part of a transistor is placed on top of a substrate, i.e., 
it has a "mesa" geometry. Transport current through 
the mesa-structure may result in large self-heating due 
to power dissipation in a small volume and a poor ther- 
mal contact to the substrate. It is well known that self- 
heating can considerably deteriorate the performance of 
high-power transistors especially at cryogenic tem- 
peratures H . Therefore, details of thermal transport in 
small mesa-structures is important for future develop- 
ment of cryoelectronics. 

High-T c superconducting (HTSC) Bi2212 mesa- 
structures arc promising candidates for cryoelectronics 
(g,|J|J^l. Due to the layered structure of Bi2212, such 
mesas represent natural stacks of atomic scale " intrinsic" 
Josephson junctions (IJJ's) Indeed, estimation of the 
stacking periodicity from Fraunhofer modulation of Fiske 
steps D yields s = 15.51, in a perfect agreement with 
the crystalline structure of Bi2212. The atomic perfec- 
tion of the naturally occurring tunnel junctions in Bi2212 
provides a reliable basis for cryo-devices giving a high de- 
gree of homogeneity and reproducibility. IJJ's with their 
record-large I c R n values ~ 5 — 10mV^ are particularly at- 
tractive for high frequency applications and have 
become a powerful tool for fundamental studies on HTSC 
materials [ ^0|JTl]Jl2]| . However, self-heating may become 
an acute problem in Bi2212-mesas |]l3| , |p4] , ^5| due to low 
thermal conductivity. 

In this letter we study self-heating in small mesa- 
structures both analytically and numerically. An ana- 
lytic solution for a three-dimensional temperature dis- 
tribution shows that self-heating is proportional to the 
length of the mesa and that it can be negligible in small 
Bi2212-mesas with a small number of IJJ's, N. 

The problem of self-heating in Bi2212 mesas has long 
been anticipated. A pronounced back-bending (negative 
differential resistance) of current-voltage characteristics 
(IVC's) is typically observed at large bias. The back- 
bending increases (see e.g., Fig. 10a in Ref. |lfj[ ) and 



the voltage separation between quasiparticle branches in 
IVC's decreases with dissipated power (e.g., see Fig. 3 
in Ref. Qj). Such a behavior is in qualitative agree- 
ment with progressive self- heating and/or nonequilib- 
rium quasiparticle injection in the mesas. A substan- 
tial self-heating was found in short-pulse experiments |TJ] 
and direct mesa temperature measurements [ p"5| . 

However, recently it was argued that self-heating is less 
important in small mesas PjTIj). It was noted that the 
back-bending of IVC's vanishes and the spacing between 
quasiparticle branches remains constant (see eg., Fig. 3 
in Ref. JllJ ) when the in-plane mesa size is decreased to a 
few /jm. Seemingly, this contradicts a simple estimation 
of the overheating 

AT ~ qt c /n c , (1) 

which suggests that the self-heating should be inde- 
pendent on the area of the mesa as the dissipated power 
density q ~ j c V is constant for mesas with the same j c . 
Here t c is the thickness and k c is the transverse thermal 
conductivity of the Bi2212 single crystal, and j c is the 
c-axis critical current density. 

To show that the self-heating does decrease in small 
structures let us consider a disk-like mesa with in-plane 
radius a, on top of a bulk single crystal, see Fig. 1. We 
assume that the stationary heat conduction is described 
by a heat diffusion equation: 

V • kVT = -V-q, (2) 

with anisotropic thermal conductivity, K a t along the 
ab-planes and k c in the transverse direction (c— axis). As- 
suming that the power is dissipated only inside the mesa, 
boundary conditions become n c dT / dz = —q at the bot- 
tom of the mesa, z — 0, \r\ < a, and T = Tq, at the 
bottom of the crystal, z = —t c . 

In the crystal, Eq. || is reduced to the Laplace equation 
by substitution z* = z-\f K a {,/ k c . Then the problem can 
be mapped to that of an electrical potential of a charged 
conductor. The exact solution is known for a 3-axis el- 
lipsoid |lq| . For a circular mesa on top of a semi- infinite 
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FIG. 1. Results of numerical calcu- 

lation with self-consistent re(T). There are 19 isotherms in 
the temperature range 10-41 K. 

single crystal t c —* oo the solution is 
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where r* 2 = r 2 + z* 2 and AT = T - T . 

A solution for a finite t c can be obtained by an image 
method: 



AT{t c ) = ATqo (r, z) - AT^r, 2t c - z). 



(4) 



From Eqs. (3,4), the self-heating of the mesa, z = 0, 
\r\ < a is 



AT m = 
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Now it is clear that Eq. [l] is only justified for a 
2t c y / K a b/ Kc ~ 60 /im for a typical t c ~ IQptm and 
Kab/Kc — 10, i.e., in none of practical cases. 

For smaller mesas, 
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which demonstrates that the self-heating decreases lin- 
early with a decrease of the in-plane size of the mesa. 

To include the temperature dependence, k(T), we 
solved self-consistently Eq. |^ using FEMLAB package 
0, taking empirical k ab {T) » (1.06 x 10~ 3 T - 7.57 x 
10" 6 T 2 ) Wcm^K- 1 @ and n ab / ' K c ~ 10 @. Results 
of numerical simulations for a disk-like mesa with the 
radius a — 5 fim on top of a single crystal 1 mm in di- 
ameter with t c = 10 fim and Tq — 10 K are shown in 
Fig.l for a semi-space, r > 0. The numerical solution 
also revealed a linear decrease of AT m with a, but the 



self-heating was smaller due to an increase of k with T. 
For example, in Fig.l the maximum temperature of the 
mesa is ~ 41 K for an empirical fc Q &(T) and ~ 67 K for 
the case of k a b = fc a b(To) = const. 

So far we assumed that the heat is sinking only into 
the pedestal through the bottom surface of the mesa. 
In reality there is always a metallic electrode on top of 
the mesa, which can help to cool the mesa. Let's as- 
sume that the electrode with the thickness t e , width w, 
length L e and thermal conductivity « e , is attached to 
the top surface of the mesa. The electrode is separated 
from the Bi2212 by an insulating layer with the thick- 
ness ti and the thermal conductivity k^. Even though 
the thermal conductance of the electrode, n e wt e /L e , is 
negligible due to the small cross-section and finite length, 
it may yet provide effective spreading of the heat power. 
Effectiveness of such a heat-spreading metallic layer on 
top of the mesa is known from the study of high-power 
semiconducting transistors 0] . To estimate the effective 
thermal resistance Rth = AT/Q, we once again use the 
electric analog. The equivalent circuit is shown in Ln- 
sset of Fig. 1. Here R m ~ (4aK a (,) _1 is the thermal re- 
sistance of the pedestal, see Eq.(6), R e = Ax/(wt e n e ) 
and Ri — ti/(wA.XKi) are thermal resistances of an el- 
ement Ax of the heat-spreading layer through the elec- 
trode and the insulating layer, respectively. Because the 
thickness of the insulating layer is small, ti ~ 0.1/xm, 
we used Eq.(l) for estimation of Ri. A straightforward 
calculation yields 
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where V R e Ri — [ti/(t e w 2 K e Ki)] is the effective 
thermal resistance through the electrode and insulator. 
Taking typical values of = O.lWcm^ 1 !^^ 1 (CaF2 at 
low T), K e = lWcm~ 1 K- 1 (Au), t e = U and w = 2a 
we obtain that \/R e Ri can be an order of magnitude less 
than R m . Even though this is a rough estimation, it 
shows that the parallel heat channel through the elec- 
trode can indeed be significant [^o). However, whatever 
heat channel is dominating, the self-heating at a given 
q decreases with the in-plane size. Such a conclusion 
is in agreement with earlier observations of better ther- 
mal properties of small multi-emitter transistors as com- 
pared to single-emitter transistors with the same total 
area g§. 

In Fig. 2, an experimental IVC at To = 4.2if is shown 
for a ~ 2 x 2[im 2 Bi2212 mesa, containing a large amount 
of IJJ's, N - 200. The Bi2212 is slightly overdoped with 
the superconducting critical temperature, T c ~ 90K. 
This data is an example of extremely large power dis- 
sipation in a small mesa. Indeed, because of large N, 
q is more than 5 • 10 4 W/cm 2 at the largest bias. The 
self-heating is seen from the back-bending of the IVC at 
large bias. On the other hand, multiple quasiparticle 



(7) 



2 




-6 -4 -2 2 4 6 



V(Volt) 

FIG. 2. Experimental IVC's of 2 x 2um 2 Bi2212 mesa with 
N ~ 200 at T = &.2K. a) IVC's of 3.5 x 7.5fim 2 mesa with 
N = 10 at T = 4.9K, 40K, 60 K, 82.3K and 121 K. b) Tem- 
perature dependence of the superconducting gap voltage. 

branches due to sequential switching of IJJ's into the re- 
sistive state are clearly visible up to the largest current, 
indicating that the mesa is not overheated above T c even 
at such high q. 

In inset a) of Fig. 2, IVC's at different T = 4.9K, 40K, 
60 K, 82. 3K and 121 K, are shown for another mesa with 
small N = 10 and about an order of magnitude lower 
q. Estimations from Eq.(7) yield that the self- heating in 
this case is 4-5 times less than for the case shown in Fig. 
2, i.e., AT should not exceed 20K at T = A.9K [§o|. 
Indeed, the IVC's are similar to those obtained by short 
pulse measurements |l4| with no sign of back-bending. 
The knee in the IVC's represents the sum-gap voltage, 
V g = 2NA/e, where A is the superconducting energy 
gap. The temperature dependence of V g per junction for 
the case of small self-heating is shown in inset b) of Fig. 
2. We may use this V g (T) dependence for calibration of 
self-heating at large q. For the case of Fig. 2, the V g per 
junction reduces to ~ 35mV at large bias, compared to 
~ 52mV in inset b). From this we estimate AT ~ 60K 
at V g ~ 7V in Fig. 2. Consequently AT at V g is - 15K 
for the IVC at T = 4.9Jf in inset a) of Fig.2. 

Our model uses an assumption that the heat escape 
from the mesa occurs via diffusion into the bulk of the 
single crystal. This assumption may become invalid when 
the mesa height h is 100 — 200 A only, because the mean 
free path of low-frequency phonons, which make an es- 
sential contribution to the heat transport at low tempera- 
tures, can be l p ~ 1 /im 3> h [ pl| . Such (non-equilibrium) 
phonons created in one of the IJJ's would shoot through 
the whole thickness of the mesa before releasing their en- 



ergy via the electron-phonon and phonon-phonon interac- 
tions already being deep in the bulk of the single crystal. 
In other words there is an additional ballistic heat chan- 
nel, which would reduce the effective heat-power density 
in the mesa. Thus our heat-diffusion estimations should 
be rather considered as an upper limit of the overheating. 

In conclusion, we have shown analytically and numer- 
ically that self-heating decreases in small mesa struc- 
tures. Therefore, reducing the mesa size is an effective 
way of avoiding self-heating. We have demonstrated that 
self-heating can be made tolerable in small HTSC mesa- 
structures even at high dissipation power. This is im- 
portant for future applications of HTSC mesas in cryo- 
electronics and fundamental studies of HTSC. In general, 
the obtained analytic solution is useful for better under- 
standing of thermal properties of small microelectronic 
elements. 
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